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L- INTRODUCTION 


Conventional ultrasonic testing CUT) techniques for the nondestructive 
evaluation (NDE) of macrostructural defect states of materials is generally 
divided into two categories; namely, the pulse- echo and through- transmission 
techniques III . In the pulse- echo technique, the transmitting and receiving 
transducers, which may be combined into a single transducer, are coupled to 
the same face of the test structure and the reflected wave field (echo) from 
the defect is analyzed in the time domain C13 or in the frequency domain C23 . 
The echo from a defect is affected by the defect type, shape, size and 
orientation, and diffraction of the incident ultrasonic wave. The received 
signal gives an accurate measure of the defect size only in a limited number 
of cases. Information on the defect has to be obtained through interpretation 
of the signals received at a variety of receiving transducer positions C3I . 

In the through- transmission technique, the transmitting and receiving 
transducers are coupled to the opposite faces of the test structure and the 
transmitted wave field is analyzed. Although the through- transmission 
technique is less subject to errors due to defect orientation and scattering, 
it still suffers from deficiencies such as defect transparancy C3] . 

For the characterization of cracks perpendicular to the faces of a test 
structure, special separate transmitting and receiving transducers with an 
inclined beam angle may be used and the principles of conventional pulse- echo 
or through- transmission techniques may be applied. These techniques are 
referred to as tandem or pitch- and- catch technique Cl, 3, 43. 

In another technique for the characterization of such cracks, the time 
delays of the diffracted waves from the edges of the crack may be used C53 . If 
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the crack is opened to the surface (surface- breaking crack) and access to 
that surface is available, Rayleigh waves may also be utilized to characterize 
the crack £3,6,7 3. 

All the above techniques are limited to the analysis of nonoverlapping 
signals. Recently Vary et al . £8,93 introduced an ultrasonic NDE parameter 
called the “stress wave factor" (SWF). Separate transmitting and receiving 
transducers are coupled to the same face of the test structure, which is 
called the SWF configuration . An input pulse having a broadband frequency 
spectrum is applied to the transmitting transducer and the number of 
oscillations (or some modifications thereof CIO]) in the output signal at the 
receiving transducer exceeding a preselected voltage threshold is defined as 
the SWF. The SWF is also valid for the analysis of the output signals with 
overlapping echoes. 

The SWF analysis has shown encouraging results in several NDE 
applications C8- 113, and its use in new applications is currently under study 
C 1 13 . 

The stress wave transmission characteristics of an isotropic elastic 
plate through which the transmitting and receiving transducers communicate in 
the SWF configuration are studied in C123 and C 1 33 . The analysis of SWF output 
signals in the frequency domain is given in £143 . From the analyses in C12- 
143, the following can be stated for the SWF signal. 

I . The SWF signal is the superposition of signals due to stress wave 
fields which multiply reflect between the top and bottom faces of the 
structure, simply reflections. Each reflection comes from a different 
direction and travels a different total ray path length. Reflections 


may be observed separately if the plate is sufficiently thick, or 
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they may overlap. 

2. More advanced signal processing techniques C141 are needed for the 
analysis of the output, especially if the reflections are overlapped. 
However, more information can be obtained by analyzing only one 
output . 

3 . The output contains information about a larger volume of the plate as 
compared with the conventional techniques, and there is flexibility 
in specifying the testing region via changing the distance between 
the transmitting and receiving transducers. 

4. The output is affected by the wave transmission properties of the 
plate in the direction parallel to its faces, the direction in which 
plates are loaded in many structures. 

5. UT by utilizing the SWF configuration does not require access to both 
sides of the test structure. 

6 . If the thickness of the structure is less than the near field 

distance C12, 131 (which depends on the transmitting transducer 
radius and the signal wavelength) the measurements in the 

conventional pulse- echo and through- transmission techniques must be 
made in the near field, or the transmitting transducer and the plate 
may have to be located apart from each other in a fluid tank. 
However, in the SWF configuartion it may be possible to get 
information in the far field although the plate is thin, because 
there may be multiple reflections which travel a total ray path 
length greater than the far field distance C131. The measurements in 
the far field are outside the interference region and, therefore, may 


be interpreted more easily. 



4 


The SWF has been generally used to characterize microstructural defect 
state of materials [8- 113. In this study the experimental SWF output signals 
generated in plates containing cracks perpendicular to its faces are analyzed 
using the signal processing techniques discussed in C143. The effects of the 
existence, lateral location and depth of the crack on the output signal 
characteristics are studied in the frequency domain. 
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2 . EXPERIMENTS 


The experimental system used in this study is characterized in C153 by 
coupling the transmitting and receiving transducers face- to- face, without a 
test specimen in between. The experimental system behaves as a linear- phase 
bandpass filter in the frequency range 0.6 to 2.3 MHz. 

A schematic of the SWF configuration showing the experimental geometry 

is given in Fig. 1. The plate has a thickness h and its sides are far from the 

transducers so that reflections from the sides do not occur during the period 

of observation . The transducers have a radius a and are located on the top 

face of the plate at a distance £ apart. The crack has a depth d which is 

perpendicular to the plate faces. It is open to the bottom face of the plate 

throughout the depth (into the plane of the sketch) of the plate. It is 

located at a distance £ from the middle axis between the transducers, where 

c 

is positive if the crack is located on the receiving transducer side of the 
middle axis as in Fig. 1. 

The plates (25 cm x 25 cm x h) used in the experiments are made of 6061- 
T6 aluminum with longitudinal (P) and shear (S) wave velocities of 0.632 and 
0.313 cm/M sec, respectively. The thicknesses of the plates and the depths of 
the cracks considered are listed in Table 1 . All the cracks have a width of 
0 . 05 cm . The codes assigned for specimens and some dimensionless quantities 
associated with the geometry of specimens and the P wavelengths 
corresponding to several frequencies f are also given in Table 1 . 

The transducers (Panametrics model V105) have a radius a of 1.125 cm, 
and the distance between the transducer axes i is 6 cm in all experiments. For 
specimens with nonzero crack depth, the location of the crack € 

c 


is changed 
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from -4 to 4 cm with in increments of 0.5 cm. 

The time history of the input signal is shown in Fig. 2. It has a -6 dB 
frequency band of 0.4 to 3 MHz. The signal in Fig. 2 and all the time domain 
signals in the subsequent figures are discrete; their sampling period being 
0.05 J*sec (sampling frequency 20 MHz). The discrete points are connected by 
linear interpolation . 

The input signal in Fig. 2 is the same in all experiments except for 
specimens with h=10 cm <C0 # Cl, C2) . For specimens with h=10 cm the input has 
the same waveform but scaled by approximately 0.3 compared with the input 
signal in Fig. 2. 
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3. PROCESSING OF SWF OUTPUT SIGNALS 


Typical SWF output signals for some cases are shown in Figs. 3-5. As in 
C13] , the reflection which is generated fay the stress wave fields which travel 
p times as a P wave and s times as an S wave through the plate while 
reflecting back and forth between the bottom and top faces of the plate is 
designated as the reflection with p,s. The time delays for some possible 
reflections are calculated using the formulae given in C123 and C 1 3D , which 
are derived by considering the multiply reflected rays constructed 
geometrically according to Snell's law; these are indicated above the output 
signals in Figs. 3, 4 and 5 with the values of p and s. 

As observed from Fig. 3, the reflections in the output signals for 
specimens AO and A1 (with h=2.4 cm) are not overlapped significantly. This 
nonoverlapping case exists because the thickness of the plate is sufficiently 
large that the difference between the time delays for the successive 
significant reflections is adequately large. On the other hand, the 
reflections in the output signals for specimens BO and B1 (with h=1.25 cm) are 
overlapped as observed from Fig. 4. Also, the output signals shown in Fig. 5 
for specimens CO and C2 (with h = 10 cm) contain nonoverlapping reflections. 

The digital signal processing techniques applied to SWF signals to 
obtain information in the frequency domain on individual reflections are 
discussed in C141 . Short- time portions of the SWF signals are considered. If 
significant reflections in a SWF signal are not overlapped, the short- time 
Fourier analysis is used. The data window, which defines the short- time SWF 
signal, includes only an individual reflection and completely excludes all the 
other reflections. Then, by calculating the fast Fourier transform (FFT) of 




the short- time SWF signal, the magnitude spectrum of the individual 
reflection is found directly. 

If reflections in a SWF signal are overlapped, the short- time 
homomorphic analysis, which is also called the short- time cepstrum analysis, 
is used. The inverse fast Fourier transform (IFFT) of the logarithmic 
magnitude spectrum (log spectrum) of a time domain signal is called the 
cepstrum. The cepstrum is again a time domain signal and is symmetric about 
zero. In order to find an estimate of the magnitude spectrum of an individual 
reflection, the cepstrum of the short- time SWF signal is multiplied by a 
cepstrum window which covers the low- time portion of the cepstrum. Then the 
FFT and exponentiation operations are applied to the windowed cepstrum, 
resulting in the estimate to the magnitude spectrum of the individual 
reflection. For the analysis to be valid, the data window and the cepstrum 
window must be specified appropriately, the conditions for which are discussed 
in CM! in detail. 

The types of the methods of analysis and the specifications of the data 
windows for finding the magnitude spectra of , some reflections with only P 
waves (s = 0) are given in Table 2 for specimens AO, BO, CO (with d=0) . The 
reflections considered for each specimen are specified by the values of p. The 
ray angles, total ray path lengths, and time delays for the reflections are 
calculated using the formulae given in C123 and C133. These are also indicated 
in Table 2. The cepstrum window used in the short- time homomorphic analysis 
is a Hamming window with an interval of -1.25 to 1.25 Msec. 

For any other specimen having the same plate thickness as any of the 
specimen AO or BO or CO the corresponding same conditions given in Table 2 are 


considered . 
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A. RESULTS AND DISCUSSIONS 

The characteristics of a log spectrum such as the peak frequency f ; 

P 

lower and upper frequencies of —3 dB frequency band f and f , respectively; 

• A 

lower and upper frequencies of -6 dB frequency band f and f , respectively; 
average resonant frequency f peak log magnitude Y p ; log magnitude at average 
resonant frequency Y^; and average log magnitude over - 3 dB frequency band 
(simply the average log magnitude) Y^ are defined in Fig. 6 on a hypothetical 
log spectrum. 

The values of the average log magnitude and the average resonant 

frequency for a particular reflection in a specimen without a crack are 

denoted by CY 3 and Cf 1 ^respectively; and these are listed in Table 3 for 
a o co 

different reflections in specimens AO, BO, and CO. The statistical behaviors 

°* C ^a^o anc * Cf c ] ° w k* c h are evaluated from three sets of repeated 

experiments, each are also given in Table 3. Notice that the frequencies 

considered in Table 1 cover the range of frequencies that appears in Table 3. 

How the existence, lateral location, and depth of the cracks in 

specimens with h=2.4 cm (AO, Al, A2, and A3) affect the average log magnitude 

the log spectra of the reflections considered is studied in Fig. 7. 

Figs. 7 (a) , 7 <b) , and 7 <c) show the behavior of Y -CY 1 for the 

a a o 

reflections with p=8, 10, and 12, respectively, versus £ for different d The 

c 

locations of the transducers and the middle axis are indicated in the figures. 

The data points at the beginning and the end of the abscissa (location of 
crack axis) in the figures correspond to the specimen without crack (AO). They 
are generated from separate experiments “under identical experimental' 1 
conditions. For the data point at the beginning of the abscissa, Y -CY ] is 
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exactly zero. For the data point at the end of the abscissa# Y^— is 
approximately zero. 

If there were no crack in the specimens with h=2.4 cm then all of the 

specimens would be identical and Y -CY 3 would be approximately zero within a 

a a o 

threshold value for all the data points in Figs. 7 (a), 7(b) and 7 (c) . The 

value of the threshold limit depends on the repeatibility of the experiments 

"under identical experimental" conditions. This value can be estimated using 

the standard deviation values for CY 3 given in Table 3. Thus# the changes in 

a o 

Y -CY 3 above the threshold limit in Figs. 7 (a), 7(b) and 7 (c) are due to 
a a o 

the the the existence of cracks. 

There are several observations which can be made from Figs. 7 (a), 7 

<b> , and 7 <c) . 

1. The change in Y -CY ] versus € is approximately symmetric about the 

a a o c 

middle axis. 

2. Y -CY 1 is negative if the crack is located between the transducer 

a a o 

axes; that is, if I * i <3 where l l denotes the "magnitude of". On 

c 

the other hand, Y -CY 3 "is generally positive" when I € I >3. 

a a o c 

3. The magnitude of the negative values of Y^-CY^l^ is increasing for 

increasing crack depth d in the region where 1 5 2. The average of 

the data points in the region with I €^1 < 2 for a particular 

reflection (p) and a particular crack depth (d) is called the average 

dB drop and it is denoted by CY 3 -CY 1 

a m a o 
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4. The average dB drop versus 100 d/h is shown in Fig. 7 <d) for the 
different reflections (with p=8, 10, and 12) considered in specimens 
with h=2.4 cm. It is observed that the average dB drop for a 
particular d/h is approximately the same for each reflection 
considered . 

Fig. 8 shows how a particular reflection (with p=4, s=0) interacts with 
a crack. It is shown in C13] that a receiving point in the far field observes 
the stress wave field for a particular reflection as a plane wave propagating 
and multiply reflecting in the directions specified by the ray which is 
constructed geometrically between the receiving point and the center of the 
transmitting transducer using Snell's law for that particular reflection. Such 
a multiply reflected ray is constructed in Fig. 8 (a). 

In order to observe the scattering effect of the crack more clearly, 
consider the hypothetical multi- layered half- space shown in Fig. 8 (b) , 
which is introduced in C13] . It is assumed that the successive back and forth 
reflections at the bottom and top faces of the original plate take place as 
transmissions through the boundaries of layers. The layers are assumed to be 
bonded together in a hypothetical way such that an incident wave in one layer 
produces no reflected waves but produces transmitted waves in the next layer 
with the transmission coefficients that are the same as the reflection 
coefficients as if the incident wave reflects at a stress- free plane 
boundary. In this description, the top layer can be thought of as the original 
layer and each subsequent layer can be thought of as the image of the layer 
immediately above it. The ray corresponding to the multiply reflected ray in 
Fig. 8 (a) is depicted in Fig. 8 (b) . Notice that the images of the crack also 


appear in the image layers in Fig. 8 (b) . 
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The description in Fig . 8 (b) then suggests that a reflection with a 

particular value of p and s "sees" the <p+s)/2 cracks which are twice the size 
of the original crack. If the cracks are located between the transducer axes 
as in Fig. 8 <b> , they direct some of the incident energy away from the 
receiving transducer, which may cause a drop in the energy of the waves 
arriving at the receiving transducer, as compared with the no- crack case. If 
the cracks are not located between or on the transducer axes, they direct some 
of the incident energy toward the receiving transducer; energy which would 
otherwise not reach the receiving transducer if there were no cracks. This may 
cause an increase in the energy of the waves arriving at the receiving 
transducer . 

As the location of a crack changes, scattered field from the crack 

changes and thus the field observed at the receiving transducer changes. For 

the crack sizes, location of transducers and frequencies considered in plates 

with h=2 .41 cm, these changes do not affect the condition that the larger the 

crack the larger the drop in Y -CY 3o observed in Figs . 7 (a) , 7 (b) , and 7 

a a 

(c) . 

The reflections with larger p <s=0> are affected by the crack with 
larger times as discussed above. On the other hand, for increasing p the angle 
between the propagation direction and crack becomes smaller, which decreases 
effect of the crack on the wave. This explains qualitatively why all the 
reflections undergo approximately the same average dB drop for a particular 
d/h, as seen in Fig. 7 (d) . 

The analysis for specimens with h=2.4 cm in Fig. 7 is repeated for 
specimens with h=1.25 cm (BO, B1 , B2, and B2) and with h=10 cm <C0, Cl, and 
C2> in Figs. 9 and 10, respectively. The data for the average dB drop versus 
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100 d/h for different groups of specimens in Figs. 7 (d) , 9(d), and 10 <d> are 
superposed in Fig. 11. The line which is fitted to the superposed data by 
using least- squares criterion [14] is also shown. 

Each data point in Fig. 11 satisfies the following condition : 

d 

100 — — = -7 {CY 3 -CY 3 3 + Error (1) 

h a m a o 

where -7 is the slope of the correlation line in Fig. 11. 

The error values in eqn, (1) for the data points in Fig. 11 are given in 

Table 4. The absolute value of the error is less than 10 in most cases, 

however it reaches a maximum value of 14.9. 

How the existence, lateral location and depth of the cracks affect the 

peak log magnitudes Y^ and the log magnitudes at average resonant frequencies 

Y are also studied. The same results shown in Figs. 7 and Figs. 9-11 are 

obtained, but they are not repeated here. 

The effect of cracks on the average resonant frequencies f is studied 

c 

for the same reflections considered for the specimens with h=2.4 cm. The 

results are shown in Fig. 12. Again, the change in f — Cf 3 versus { is 

c c o c 

symmetric about the middle axis. Here the change in f -[f 3 is larger if the 

C CO 

crack is located near (either) one of the transducer axes. However, there is 

no simple interpretation in terms of correlating the change in f -Cf 3 to 

c co 

the crack depths. 

Similar results are obtained for other groups of specimens with h=1.25 

cm and h=10 cm. Again no simple correlation between f -Cf 3 and the crack 

c c o 

depths are observed. 

Similar results are also obtained for the peak frequencies f ; lower and 

P 
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upper frequencies of -3 dB frequency band ^ and f respectively; and lower 
and upper frequencies of -6 dB frequency band f and f , respectively. Again 

W 1 

no simple correlation to the crack depths can be stated at this point . The 
effect of a crack on all these characteristics are under current 
investigation. 
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5. CONCLUSIONS AND RECOMMENDATIONS 

The stress wave factor (SWF) signal, which is the output of an 
ultrasonic testing (UT) system where the transmitting and receiving 
transducers are coupled to the same face of the test structure, was analyzed 
for the nondestructive evaluation (NDE) of plates containing perpendicular 
edge cracks. The transducers were located at the top face of the plate. The 
crack was open to the bottom face throughout the depth of the plate. 

The effects of the existence, lateral location and depth of the crack on 
the log spectra of individual reflections in the SWF signal were studied. If 
the reflections in the SWF signal were not overlapped, the short- time Fourier 
analysis was applied to obtain the log spectra of the individual reflections. 

If the reflections were overlapped, the short- time homomorphic analysis 
(cepstrum analysis) was applied. 

Aluminum plates with thicknesses h=1.25, 2.4, and 10 cm were used. 
Cracks with 100 d/h= 25, 50, and 75 in plates with h=1.25 and 2.4 cm, where d 
was the crack depth, and cracks with 100 d/h= 6 and 12 in plates with h=10 cm 
were considered . The transducers which had a radius of 1 . 125 cm were located 6 
cm apart . 

Plates with h=1.25 cm contained overlapping reflections and plates with 
h=2.4 and 10 cm contained nonoverlapping reflections. The log spectra of 
several reflections which had average resonant frequencies approximately at 
0.9, 1.3, and 1.7 MHz were analyzed. It was observed that the average log 
magnitudes of the log spectra decreased more with increasing d/h if the crack 
was located between the transducers . The average decrease in average log 


magnitudes for different crack locations between the transducers was called 
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the average dB drop. 

It was observed that the average dB drop versus 100 d/h was 
approximately linear. Also, the linearity was approximately the same for each 
reflection and each group of plates with the different thicknesses considered. 

By using the least- squares criterion, the slope of the correlation line for 
all the data was found to be 1 dB per 0.07 d/h. 

Changes in the average resonant frequencies of the log spectra were also 
observed due to the existence of a crack and to changes in the location of the 
crack. However, no simple correlation between the average resonant frequency 
and crack depth could be stated at this point. 

This study provides a quantitative analysis of the effects of the 
existence, lateral location and depth of a perpendicular crack on the SWF 
signal in plates. The analysis should be extended for cracks with different 
orientations and shapes and for different frequency ranges. Scattering 
theories should also be utilized to develop more quantitative analyses. 
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TABLE 1 Geometrical Properties of Specimens and Some Associated 
Dimansionless Quantities . 
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TABLE 2 Methods of Analysis and Specifications of Data Windows to Find 
Magnitude Spectra of Reflections With Only F Waves; and Rav Angies, 
Total Ray Path Lengths , and Time Delays of Reflections in Some 
Specimens . 


Soecimen 

Code 

Reflection 
(Value of p) 

Data Window 

Ray 

Angle 

(Degree) 

Total Ray 
Path Length 
(cm) 

Time 
Delay 
()* sec) 


Interval 0*sec> 

From 

To 


3 


31.8 

35.3 

17.35 

20.12 

31.83 

AO 

10 

STF* 

39.1 

42.6 

14.04 

24.74 

39.14 


12 


46.5 

50 

11.77 

29 42 

46.55 


12 


16.95 

34.95 

21.80 

16.16 

25.56 

BO 

18 

STH b 

28.45 

46.45 

14.93 

23.29 

36.85 


22 


35.95 

53.95 

12.31 

28.15 

44.54 


2 


31.5 

39 

16.70 

20.88 

33.04 

CO 

4 

STF* 

62 

69.5 

8.53 

40.45 

64.00 


6 


93 

100 .,5 

5.71 

60.30 

95.41 


a STF : Short-time Fourier analysis (512-point FFT is used.) 
b STH : Short-time homomorphic anaivsis <512-point FFT and IFFT are used 
in the process. The cepstrum window used in the process is a 
Hamming window having a duration from -1.25 to 1.25 **sec.) 



































TABLE 3 Values of f and M (f ) for Some Reflections in Specimens 
o o o 

Without Crack and their Statistical Behaviors. 






Mean Value 

Standard 

Mean Value 

Standard 

Specimen 

Reflection 

t 

0 

M (f ) 
o o 

of f 

o 

Deviation 

of M (f ) 
o o 

Deviation 


(Value of p> 

<MHr> 

<dB> 

(MHz) 

in f 

0 

(dB) 

in M (f ) 
o o 






(MHz) 


(dB) 


8 

1 

8 32 

0.95 

0 034 

8.54 

0.18 

AO 

10 

l . 19 

9.33 

1 16 

0.018 

9.87 

0.39 


12 

1 49 

7.38 

1.47 

0.015 

8 23 

0.62 


12 

0.97 

-1 13 

0.96 

0.010 

-1.24 

0.85 

BO 

18 

1.38 

-2.14 

1.37 

0 014 

-2 03 

0.16 


22 

1 41 

-4.38 

1.40 

0 018 

-4.38 

0.35 


2 

1 O.BB 

6 64 

0.86 

0.014 

6 68 

0.04 

CO 

4 

1.37 

16 09 

1.36 

0.006 

16 07 

0.08 


6 

1.68 

14 01 

1.65 

0.021 

14.33 

0.26 
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Fit* - 2 Tim* history of input signal. 
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Fig. 3 Time histories of output signals tor specimens with < c *° toT 


<a> AO <h-2.4 cm, d»0>, and <b> A1 (h-2.4 cm,d»0.6 cm). 
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Fig. 4 Tim® histories of output signals for specimens with f^aO for 


(a) BO (h> 1.25 cm, d=0>, and (b> B1 (ha 1.25 cm, da0.31 cm) 
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Fig. 5 Tim* histories of output signals for specimens with t for (a) CO 

c 

(h 1 0 cm, daO), and <b) C2 (h»10 cm, d»1.2 cm). 
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Fig. 7 M(f )/M <f ) for reflections with (a) p»8, (b) p«10, and (c) p*12 
ooo 

versus l for different d; and <d> 20 LoaCM(f )/M (f >3 versus d/h for 
c ooo 

different p in specimens with h*2.4 cm (AO, A1 , A2, and A3). 
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Fig. 8 <a> Plate with crack and the rav for the reflection with o=4,s 

and (b) hypothetical multi-layered medium description . 
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Fig. 10 M<f >/M (f > for reflections with <e> p«2, (b> p*4, end (c) 

^ o o o 

versus t Q for different d; end <d> 20 LogCM(f o )/M o (f o )3 versus d/h 
for different p in specimens with h»10 cm (CO, Cl, end C2>. 
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Fig. 11 Superposition of the data in Figs. 7 <d> f 9 <d> , and 10 (d) 
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line fitted to superposed data by using least-squares criterion. 
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